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Abstract 

Transformation optics constructions have allowed the design of 
cloaking devices that steer electromagnetic, acoustic and quantum 
waves around a region without penetrating it, so that this region is 
hidden from external observations. The material parameters used to 
describe these devices are anisotropic, and singular at the interface 
between the cloaked and uncloaked regions, making physical realiza- 
tion a challenge. These singular material parameters correspond to 
singular coefficient functions in the partial differential equations mod- 
eling these constructions and the presence of these singularities causes 
various mathematical problems and physical effects on the interface 
surface. 

In this paper, we analyze the two dimensional cloaking for Helmholtz 
equation when there are sources or sinks present inside the cloaked re- 
gion. In particular, we consider nonsingular approximate invisibility 
cloaks based on the truncation of the singular transformations. Using 
such truncation we analyze the limit when the approximate cloaking 
approaches the ideal cloaking. We show that, surprisingly, a non-local 
boundary condition appears on the inner cloak interface. This effect 
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in the two dimensional (or cylindrical) invisibility cloaks, which seems 
to be caused by the infinite phase velocity near the interface between 
the cloaked and uncloaked regions, is very different from the earlier 
studied behavior of the solutions in the three dimensional cloaks. 

1 Introduction 

There has recently been much activity concerning cloaking, or rendering ob- 
jects invisible to detection by electromagnetic, acoustic, or other type of 
waves or physical fields. Many suggestions to implement cloaking has been 
based on transformation optics, that is, designs of electromagnetic or acous- 
tic devices with customized effects on wave propagation, made possible by 
taking advantage of the transformation rules for the material properties of 
optics. All perfect cloaking devices based on transformation optics require 
anisotropic and singular material parameters, whether the conductivity (elec- 
trostatic) [16, 17], index of refraction (Helmholtz) [22], [10], permittivity and 
permeability (Maxwell) [29], [10], mass tensor (acoustic) [10], [6], [9], or ef- 
fective mass (Schrodinger) [13, 14, 35]. By singular material parameters, 
we mean that at least one of the eigenvalues or the values of the functions 
describing the material properties goes to zero or infinity at some points 
when the material parameters are represented in Euclidean coordinates, typ- 
ically on the interface between the cloaked and uncloaked regions. Both the 
anisotropy and singularity present serious challenges in trying to physically 
realize such theoretical plans using metamaterials. Analogous difficulties are 
encountered in the study of invisibility cloaks base on ray-theory [23] and 
plasmonic resonances [1, 26]. 

To justify the invisibility cloaking constructions, one needs to study phys- 
ically meaningful solutions of the resulting partial differential equations on 
the whole domain, including the region where material parameters become 
singular. In [10], the finite energy solutions are defined to be at least measur- 
able functions with finite energy in (degenerate) singular weighted Sobolev 
spaces, and satisfy the equations in distributional sense. 

Due to the presence of singular material parameters, or mathematically 
speaking, partial differential equations with singular coefficient functions, the 
question how the waves behave in cloaking devices near the surface where 
the material parameters are singular is complicated. Indeed, very different 
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kind of behaviors of solutions are possible: In the three dimensional case, 
it is proved in [10] that the transformation optics construction based on 
a blow up map allows cloaking with respect to time-harmonic solutions of 
the Helmholtz equation or Maxwell's equations as long as the object being 
cloaked is passive. In fact, for the Helmholtz equation, the object can be 
an active source or sink. Moreover, in [10] it is shown that the finite energy 
solutions for the Helmholtz equation in the three dimensional case satisfy 
a hidden boundary condition, namely waves inside the cloaked region sat- 
isfy the Neumann boundary conditions. For Maxwell's equations, the finite 
energy solutions inside the cloaked region need to have vanishing Cauchy 
data i.e., the hidden boundary conditions are over-determined. This leads 
to non-existence of finite solutions for Maxwell's equations with generic in- 
ternal currents [10]. Physically, this non-existence results is related to the 
so-called extraordinary boundary effects on the interface between the cloaked 
and uncloaked regions [36]. 

Another point of view in dealing with the singular anisotropic design for 
cloaking devices is to approximate the ideal cloaking parameters by nonsin- 
gular, or even nonsingular and isotropic, parameters [12, 13, 14, 18, 19, 25], 
which has its advantages in practical fabrication. In the truncation based 
nonsingular approximate cloaking for three dimensional Helmholtz equation 
[14], when it approaches the ideal cloaking, one can obtain above Neumann 
hidden boundary condition for the finite energy solution. In the nonsingular 
and isotropic approximate cloaking, one can obtain different types of Robin 
boundary conditions by varying slightly the way how the approximative cloak 
in constructed, see [12, 13, 14]. 

Similarly, for Maxwell's equations it has been studied how the approximate 
cloak behave on the limit when the approximate cloaks approach the ideal 
one [25]. We note that for Maxwell's equations there are various suggestions 
what kind of limiting cloaks are possible in three dimensions. These sug- 
gestions are based on constructions where additional layers (e.g. perfectly 
conduction layer) is attached inside the cloak [10] or where the ideal cloak 
corresponds to some of the possible self-adjoint extensions of Maxwell's equa- 
tions [33, 34]. In the two dimensional or cylindrical cloaking construction for 
Maxwell's equations, the eigenvalues of permittivity and the permeability of 
the cloaking medium do not only contain eigenvalues approaching to zero 
(as in 3D) but also some of the eigenvalues approach infinity at the cloaking 
interface. Then, the electric flux density D and magnetic flux density B may 
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blow up even when there are no sources inside the cloak and an incident plane 
wave is scatters from the cloak, see [11]. However, if a soft-hard (SH)-surface 
is included inside the cloak, the solutions behave well. 

These above examples show how different behavior the solutions may have, 
in different type of cloaking devices, near the interface between the cloaked 
and uncloaked regions. 

In this paper, we analyze the two dimensional cloaking for Helmholtz equa- 
tion when there is a point source inside the cloaked region. We start with 
the nonsingular approximate cloaking based on the truncation of the singu- 
lar transformation. Taking the limit when the approximate cloaks approach 
the ideal cloak, we show that a non-local boundary condition appears on the 
inner cloak interface. This type of boundary behavior is very different from 
that the solutions have in three dimensional case discussed in [14, 36]. The 
main result is formulated as Theorem 3.2. Physically speaking, such non- 
local boundary condition is possible due to the fact that the phase velocity 
of the waves in the invisibility cloak approaches infinity near the interface 
between the cloaked and uncloaked regions, even though the group veloc- 
ity stays finite, see [7]. We note that as the most important experimental 
implementations of invisibility cloaks [30] have been based on cylindrical 
cloaks, the appearance of such boundary condition could also be studied in 
the present experimental configurations, at least on micro-wave frequencies. 
We also study the eigenvalues, i.e., resonances inside the ideal cloak corre- 
sponding to the non-local boundary condition. As these eigenvalues play an 
essential role in the study of almost trapped states [13] and in the develop- 
ment of the invisible sensors [2, 15], such resonances can be used to study 
analogous constructions in the cylindrical geometry. 

The rest of the paper is organized as following. In Section 2, we consider 
the basics on transformation optics in the electrostatic setting and the ideal 
acoustic cloak for the two dimensional Helmholtz equation. Section 3 is 
devoted to the nonsingular approximate acoustic cloaking construction and 
analysis of behaviors of acoustic waves as it approaches the ideal cloaking. 
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2 Perfect acoustic cloaking 



2.1 Background: electrostatics 

Our analysis is closely related to the inverse problem for electrostatics, or 
Calderon's conductivity problem [3, 5, 27, 28, 32]. Let Q C M. d be a domain, 
at the boundary of which electrostatic measurements are to be made, and 
denote by cr(x) the anisotropic conductivity within. In the absence of sources, 
an electrostatic potential u satisfies a divergence form equation, 

V-(jVu = (1) 

on Q. To uniquely fix the solution u it is enough to give its value, /, on the 
boundary. In the idealized case, one measures, for all voltage distributions 
u\an = f on the boundary the corresponding current fluxes, u-aVu, where 
v is the exterior unit normal to dfl. Mathematically this amounts to the 
knowledge of the Dirichlet-Neumann (DN) map, A CT . corresponding to cr, 
i.e., the map taking the Dirichlet boundary values of the solution to (1) to 
the corresponding Neumann boundary values, 

K ■ u\ d n h> u-aVu\ dQ . 

If F : Q, — > fl, F = (F 1 , . . . , F d ), is a diffeomorphism with F\gn 
then by making the change of variables y = F(x) and setting u 
we obtain 

V ■ aVv = 0, 
where a = F*a is the push forward of a in F, 

This can be used to show that 

Thus, there is a large (infinite-dimensional) family of conductivities which 
all give rise to the same electrostatic measurements at the boundary. This 
observation is due to Luc Tartar (see [20] for an account.) Calderon's inverse 
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problem for anisotropic conductivities is then the question of whether two 
conductivities with the same DN operator must be push-forwards of each 
other. There are a number of positive results in this direction in two dimen- 
sions [4, 21, 24, 31] , but it was shown in [16, 17] in three dimensions and in 
[19] two dimensions that, if one allows singular maps, then in fact there are 
counterexamples, i.e., conductivities that are undetectable to electrostatic 
measurements at the boundary. 

From now on, we will restrict ourselves to the two dimensional case. For 
each R > 0, let Br = {\x\ < R} and E# = {\x\ = R} be the central ball 
and sphere of radius R, resp., in M 2 , and let O = (0, 0, 0) denote the origin. 
To construct an invisibility cloak, for simplicity we use the specific singular 
coordinate transformation F : M 2 \{(3} — > M 2 \i? 1 , given by 

f V, for \y\ > 2, 

* = *W:=| + £orO<M<2. < 4 > 

Letting ao = 1 be the homogeneous isotropic conductivity on I 2 , F then 
defines a conductivity a on M 2 \5! by the formula 

a* k (x) := (F*a y k (x), (5) 

cf. (3). More explicitly, the matrix a = [cr J,fc ] 2 fc=1 is of the form 

a(x) = Jiz^n(x) + — pi— (J-n(x)), K|x|<2, 



where U(x) : 1R 3 — > M 3 is the projection to the radial direction, defined by 

U(x)v = (v~) A, (6) 

V \ X \J fI 

i.e., n(x) is represented by the matrix {x^xx*, cf. [19]. 

One sees that <j(x) is singular at Ei, that is, in interface between the cloaked 
and uncloaked regions, as one of its eigenvalues, namely the one correspond- 
ing to the radial direction, tends to 0, and the other tends to oo as \x\ 
approaching 1 + . We extend a to B\ as an arbitrary smooth, nonsingu- 
lar (bounded from above and below) conductivity there. Let Q = the 
conductivity a is then a cloaking conductivity on Q, as it is indistinguish- 
able from cr , vis-a-vis electrostatic boundary measurements of electrostatic 
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potentials (treated rigorously as bounded, distributional solutions of the de- 
generate elliptic boundary value problem corresponding to a [16, 17]). 

A transformation optics construction based on the above blow up map F 
was proposed in Pendry, Schurig and Smith [29] to cloak the region B\ in 
IR 3 from observation by electromagnetic waves at a positive frequency; see 
also Leonhardt [22] for a related proposition for the Helmholtz equation in 
IR 2 based on the use of several leaves of an Riemannian sufrace. 



2.2 Ideal cloaking for the Helmholtz equation with in- 
terior sources 

We consider the Helmholtz equation, with source term p, of the form 

AV • ffVw + co 2 u = p(x) on f2 (7) 

corresponding cloaking medium with the mass tensor and the bulk modulus 
given by 



^■jk 



where (a , A ) corresponds to homogeneous vacuum space and (cr a , A a ) are 
arbitrary smooth, nondegenerate medium in cloaked region B\. The push- 
forward of tensor F*o"o is defined by (3) and -F*Ao by 

(F*\ )(x) := [det(DF)X ]oF- 1 (x). 

More specifically, if (cxo, Ao) = (/, 1), one has for 1 < \x\ < 2, 

a {x) = M^n(x) + 1 -pL(/-n(z)) J X(x) 
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are both singular at the cloaking surface S := Ei = dB\. 

In the next section, we study the regularized approximate cloaking and obtain 
the behavior of acoustic waves in a singular medium by taking the limit of 
waves propagating in the nonsingular medium. 
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3 Nonsingular approximate cloaking for the 
Helmholtz equation with interior sources 



At the present time, we assume that a and A be homogeneous isotropic in 
side £>!, i.e., (a, A) = (a a 5^ k , A a ) with a a and A a arbitrary positive constants. 

To start, let 1 < R < 2, p = 2(R— 1) and introduce the coordinate transfor- 
mation F R :R 2 \B p ->R 2 \B R , 

f V, for \y\ > 2, 

X : = F ^> = { (l + M) ft, for p < \y\ < 2, 

cf. [12, 13, 14, 18, 19, 25]. We define the corresponding approximate mass 
tensor a R and bulk modulus \ R as 



a^ k (x) for \x\ > R, \ / \ _ J M x ) f° r M > ^> / Q \ 
a> fe for |x| < i?. R ^ ~ \ A a for < i?, ^ 

where a- 7 ^ and A are as in (8). Note that then a^ k {x) = ((Fr) ^ aoY* (x) 
and A(x) = ((Fr)^X ) (x) for \x\ > R. Observe that, for each R > 1, the 
medium is nonsingular, i.e., is bounded from above and below with, however, 
the lower bound going to 0, and the upper bound going to oo as R — > 1 + . 
Consider the solutions of 

(AflV- ctrV + uj 2 )u r : =p in Q 
u R \an —f, 

As a R and X R are now non-singular everywhere on Q, we have the standard 
transmission conditions on H R ■ := {x : \x\ — R}, 

Ur\x r + = u R \x R -, 
e r - (TrVur\x r + = e r - a R Vu R \x, R -, 

where e r is the radial unit vector and ± indicates when the trace on H R is 
computed as the limit r — y R^ 1 . 

Let Q, = B3. Then u R defines two functions v R such that 



u R (x) 



v R (F R \x)), for R < \x\ < 3, 
v R (x), for \x\ < R, 
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and v]i satisfy 

(V 2 + oj 2 )v+(y) = p(F R (y)) in p < \y\ < 3, 

v R\dB 3 = f, 

and 



(12) 



(V 2 + k 2 uj 2 )v^(x) = K 2 p(x), in \x\ < R. (13) 

where k 2 = (cr a A a ) _1 is a constant. Moreover, if we assume u 2 is not an 
eigenvalue of the transmission problem, then by the transformation law we 
have 

e r -a R Vu R \ m = e r -Vv%\ m . 

This implies that the DN-map K aR) \ R at dQ for the approximate cloaking 
medium (9) is the same as the DN-map at dfl, denoted by A p , of a nearly 
vacuum domain with a small inclusion present in B p . 

Next, using polar coordinates (r,6), r = \y\, and (r,6), r — \x\, the trans- 
mission conditions (11) on the surface Y, R yield 

pd r Vft(p,8) = kR drV^R, 9) . 



(14) 



For simplicity, we analyze the case of when the source is supported at origin 
inside the cloak and has the form 

p(x) = n~ 2 ^ q a d"5 {x), (15) 

\a\<N 

where 5q is the Dirac delta function at origin and q a G C, i.e., there is a 
(possibly quite strong) point source in the cloaked region. The Helmholtz 
equation (13) on the entire space M 2 , with the above point source and the 
standard radiation condition, would give rise to a radiating wave w G V(M. 2 ), 
w\ R * m G C°°(R 2 \ {0}), given by 

N 

(1)/ ,ZT\„in9 



w{r, U) 

n=-N 



J2 Pn H^(K^e m9 . (16) 



where H^(z) and J\ n \(z) denote the Hankel and Bessel functions, see [8], 
and p n G C. The analysis below can be generalized for more general sources, 
but this generalization will be considered elsewhere. 
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In Br\Br the function v R (x) differs from w by a solution to the homoge- 
neous equation of (13), and thus for r G (Ro, R) 

N 
n=-N 

with yet undefined a n = a n (K,cu; R). Similarly, for p < r < 3, 

N 



Vn(.r, 9) = (^|n|V) + b n J\ nl (ur))e md 

n=~N 



with as yet unspecified b n = b n (n,u; R) and c n = c n (n,u; R). 
Rewriting the boundary value / on dQ as 

N 



f(0) = £ U 



ind 



n=—N 



we obtain, together with transmission conditions (14), the following equations 



for a n , b n and c n : 



fn = b n J H (3cu) + c n H^(3u), (17) 
a n J\ n \{KU)R) +p n H^(KuR) = b n J\ n \{up) + c n H^(up), (18) 



KR(Kua n (J\ n \)' (kuR) + Kujp n (H^)' (kujR)) 
= p(b n u(J\ n \)'(ujp) + wc^if^/^p)). 



(19) 



Solve for a n and c n from (18)-(19) in terms of p n and 6„, and use the solutions 
obtained and the equation (17) to solve for b n in terms of f n and p n . This 
yields 

~ -t(1)/ 



■(/„ + s n H\ n (3uj)p n ) 



J\ n \(3u) + s n HW3u) 
_ (20) 

Cn S n n S n p n , 
(^n tnb-n t n p n 
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where 



s n = — {pJ| n |(reo;i?) J' n \(u)p) - n 2 RJ[ n]i (KujR)J\ n \(ujp)} 



tn = ~5~ {p H \n\( UJ P) J \n\( UJ P) ~ P H \n\ ( U P) J H P) } , 

s n = Tr {k 2 RH^'{kuR)J h {kcoR) - K 2 RJ' lnl (KuR)H^(KcuR)} , 

with D n the common denominator given by 

D n = K 2 RJ' lnl ( K uR)H^(up) - P J h {kujR){H^)'{uj P ). 

3.1 Resonances inside the cloak 

Suppose that the boundary data vanishes, i.e., / = 0. Then by (20), we have 

r(i) 

>l 



J H {3u) + s n H^(3u) 



Vr, 



Therefore, one can show 



(t n S n ~ t n 8 n )H^f(3u) - t n J\ n \(3u) 

J H (3u) + s n H$(3u) 



I 

Vn ■= ~^Vn (21) 



k 2 RH™ (KuR)h - pH™ {kuR)1 2 a 



pJ\ n \(n(jjR)l2 — K 2 RJ! n ,(KwR)li B. 



where 



h = J\ nl {up)H$(3u) - H^{up)J H {?» 



\n\ "|n| 



h = J' H (up)H^(3u) - H^\up)J H {3u) 



(22) 



11 



Since for small arguments i < 1, 



1 /x\ n 

J|n| ^ ~ \2' n - 1, J ' H ^ 

2i 



- 2 X 



n = 0, 



x n 1 n > 1, 



2™(n-l)! 



fln(x/2) n = 0, 

_*frt=iil2 n x- n n > 1, |n| W ~ 1 in\2 n x- n ~ l /-K n > 1 



2ix- 1 /n 



n = 0, 
1, 

(23) 



when n > 1, where we denote / ~ g if / — g = o(g) as x — )■ 0, and 

i2 n w -"-l( n _ 1)1 



7T 

-i2 n w- n - 1 (n - 1)! 



J, n , (3w) ujk 2 RH\ 1) '(kujR) + nif^ (kwA) 



n| 



7T 



J| n |(3w) [un 2 R J( n < (kojR) + nJ\ n \{KUjR)\ p n . 



Now we observe that \a n \ — > oo as R — > 1 + if 

[wft 2 ii( J| n |)'(«;a;i2) + nJ\ n \(KuR)] 
Note that if (25) holds, then 



0. 



(24) 



(25) 



R=l 



ujk z R(h\^)\kujR) + nH^HnuR) 



(26) 



This implies that if u; is such that (25) is satisfied by functions J\ n \ and HP] 
for some n, then there are sources p for which for any with < R\ < R2 < 1 
the H l (B R2 \ 5 Kl )-norm of the solution restricted to i?ij 2 \ B Rl goes to 
infinity (i.e., a resonance happens in the cloaked region) as R — > 1 + (i.e., 
p — > 0). We remark that condition (25) is equivalent to that the function 
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V±n(r,0) := J| n |(/sxjr)e 
satisfies the boundary value problem 

(A + kV)V = 



(27) 



in Bi, 



[nrd r V + (-d 2 ) 1/2 V] 



0. 



(28) 



r=lH 
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Roughly speaking, this means that the ideal cloak has resonance solutions 
inside the cloaked region which satisfy a non-local boundary condition. 



Next we consider the frequencies cu for which 



\ojk 2 R{ J\ n \)' (kojR) + nJ\ n \(nujR 



^0, 



R=l 



for all nSZ. 



(29) 



J\n\(ur) 



r=3 



7^0, 



3.2 Non-local boundary condition with non-resonant 
frequencies 

In the following, we show that when we have in B-2\B^ 1 < R < 2 the 
approximative cloaking material parameters, then for the non-resonant fre- 
quencies u, the boundary condition in (28) holds as well as R — > 1 + . We 
start with the following lemma. 

Lemma 3.1. Assume that in Q = B 3 we have the material parameters 
((Jji,\r). Moreover, suppose u is such that (29). When R > 1 is suffi- 
ciently close to 1, then for any source p of the form (15) with N > and 
for f G H 1 ^ 2 (d^l) the Helmholtz equation (10) has a unique solution ur. 
Moreover, as R — > 1 + , Fourier coefficients of the solution v R := ur\b r in the 
cloaked region, 



v R,n( r ) 



2- 



-ind 



v R {r,9)d9 



satisfy 



for all \n\ < N. 
Proof. We have 



\im(Krd r v Rn (r) + n 2 v Rn (r)) 



r=R 



(30) 



N 



^ {a n J\ n \(nujr) + p n H^(nur))e 
A 



in9 



n=-N 
N 

£ 

n=-N 



B, 



Kur) 



inO 
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where as in (21) 

A n = K 2 RH^'{KUjR)h - P H^(kcoR)1 2 
B n = pJ\ n \(nu)R)l2 — k 2 RJL\(kujR)Ii. 

Denote 



$ ™( r ) := ^T J \n\( KUJr ) + H \n\( KUJr )- 

Apparently $(r, 9) := $ n (r)e m£l satisfies the Helmholtz equation 

(A + = in B R . 

To check (30), it is sufficient to show 

rKd r $ n {r) + n$ n (r)\ R as R -> 1 + . 

Indeed, by 

<9 r $ n (r) = -^kujJLAkut) + (kw), 
and _B„ 7^ 0, we have 
m<9 r $ re (r) + n$ n (r)| r=i? 



r K + o;pZ 2 ) («wi2)J w (/ca;i2) - J' h {kujR)H\^{kuR) 

tin * 

where /i and Z 2 are given by (22). 
For n > 1, as p — > + , 

nZi + wp/ 2 =^P f J| fl |_i(w/9)fr| ( J Jj , (3c*;) - ^^(w/o) Jj«|(3a;) 
=0{p~ n+2 ). 

Combining (32), (33) and (24), one has 

rK(9 r $ n (r) +n$ n (r)| r=fl = 0(p 2 ) as i? -> 1 + (i.e. p -> 0" 

which proves (31). 
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For n — 0, from (23), one has 

A = - ^-^H^'(kcuR)J (3cu) In (^f) + ^RH^'{kwR)H^{3u) 
ix V 2 / 

+ —H^(kuR)J (3u) + O(p), 

1TUJ 

B =^^f (KuR)J (3u)]n (^P\ - k 2 RJ^kuR)H^(3u) 

7T V 2 / 



Therefore, 



2i 

— J (kuR)J (3u) + 0{p). 

7TU 



d r $ (R) = ^kloJq(kloR) + kujH^' (kujR) 



has denominator Bq and numerator 

koj[A Jo(kuR) + B Hjp' (kuR)] = —W n {nujR) + O(p), 

7T 

where W n (a;) = Hq 1 \x)J[ ) (x) — (x)Jq(x) This implies 

V ; kRJ^kuR) In (!f ) P 
i.e., the boundary condition (31) is satisfied for n = and moreover, 

rK0 r $ o (r)| r=fl = o(— ^ ) as R 1 + . (35) 



□ 

Now we are ready to prove our main result for the limit of the waves ur cor- 
responding to the approximate cloaking medium as R — > 1 + . The obtained 
limit can be considered as a model for solutions in an ideal (i.e. perfect) 
cylindrical invisibility cloak. 

Theorem 3.2. Let u be such that (29) is satisfied. Suppose ur is the solution 
of (10) with f = and source p is of the form (15) with N > Then as 
R — > 1 + , ur converges almost everywhere in B% to the limit u\ satisfying 



(A + k u )u\ = u p in B±, 
Krd T u x + {-dl) xl2 u 1 \ QBi = 0, 



(36) 
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and 



m i|b 2 \^ = °- (37) 



Proof. Let < Ro < R\ < R. Using Lemma 3.1 we see that the wave u R and 
its derivatives of any order converge uniformly in compact subsets of Bi\B Ro 
to 



N 
n=-N 



1 = a nJ\n\(Ku7)e ine + 



as R — > 1. Moreover, we see U\ G C°°(Bi \ B Rl ) and that U\ satisfies (38). 
As u R satisfy Dirichlet problems 

(V 2 + k 2 u 2 )u r ; = K 2 p in B Rl , 
u r \ 9Bri =h R , { ) 

where h R converge in C k (dB Rl ) to ui\q Rr for any k as R — > 1, we see that 
u R converges pointwisely to U\ in B Rl \ {0}. 

We also see that u R converge pointwisely to zero in B 3 \ B\ using the fact 
that 

c n (p) = 0(p n+l ), b n (p) = 0(p n+1 ) as p -> + . 
This proves that (37) holds. □ 

We conclude our discussion by remarking that our analysis also explains 
the limit of approximate electromagnetic cloaks in the cylindrical case with 
TE/TM polarized incoming waves, as the solutions of Maxwell's equations 
in this case satisfy the Helmholtz equation. 
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Two dimensional invisibility cloaking for 
Helmholtz equation and non-local boundary 

conditions 

Matti Lassas* and Ting Zhou^ 



Abstract 

Transformation optics constructions have allowed the design of 
cloaking devices that steer electromagnetic, acoustic and quantum 
parameters waves around a region without penetrating it, so that 
this region is hidden from external observations. The material pa- 
rameters used to describe these devices are anisotropic, and singular 
at the interface between the cloaked and uncloaked regions, making 
physical realization a challenge. These singular material parameters 
correspond to singular coefficient functions in the partial differential 
equations modeling these constructions and the presence of these sin- 
gularities causes various mathematical problems and physical effects 
on the interface surface. 

In this paper, we analyze the two dimensional cloaking for Helmholtz 
equation when there are sources or sinks present inside the cloaked re- 
gion. In particular, we consider nonsingular approximate invisibility 
cloaks based on the truncation of the singular transformations. Using 
such truncation we analyze the limit when the approximate cloaking 
approaches the ideal cloaking. We show that, surprisingly, a non-local 
boundary condition appears on the inner cloak interface. This effect in 
the two dimensional (or cylindrical) invisibility cloaks, which seems to 
be caused by the infinite phase velocity near the interface between the 
cloaked and uncloaked regions, is very different to the earlier studied 
behavior of the solutions in the three dimensional cloaks. 

* Institute of Mathematics, Helsinki University of Technology, FIN-02015, Finland. 
^Department of Mathematics, University of Washington, Seattle, WA 98195, US. 
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1 Introduction 



There has recently been much activity concerning cloaking, or rendering ob- 
jects invisible to detection by electromagnetic, acoustic, or other type of 
waves or physical fields. Many suggestions to implement cloaking has been 
based on transformation optics, that is, designs of electromagnetic or acous- 
tic devices with customized effects on wave propagation, made possible by 
taking advantage of the transformation rules for the material properties of 
optics. All perfect cloaking devices based on transformation optics require 
anisotropic and singular material parameters, whether the conductivity (elec- 
trostatic) [16, 17], index of refraction (Helmholtz) [23], [10], permittivity and 
permeability (Maxwell) [32], [10], mass tensor (acoustic) [10], [6], [9], or ef- 
fective mass (Schrodinger) [13, 14, 38]. By singular material parameters, 
we mean that at least one of the eigenvalues or the values of the functions 
describing the material properties goes to zero or infinity at some points 
when the material parameters are represented in Euclidean coordinates, typ- 
ically on the interface between the cloaked and uncloaked regions. Both the 
anisotropy and singularity present serious challenges in trying to physically 
realize such theoretical plans using metamaterials. Analogous difficulties are 
encountered in the study of invisibility cloaks base on ray-theory [24] and 
plasmonic resonances [1, 28]. 

To justify the invisibility cloaking constructions, one needs to study phys- 
ically meaningful solutions of the resulting partial differential equations on 
the whole domain, including the region where material parameters become 
singular. In [10], the finite energy solutions are defined to be at least measur- 
able functions with finite energy in (degenerate) singular weighted Sobolev 
spaces, and satisfy the equations in distributional sense. 

Due to the presence of singular material parameters, or mathematically 
speaking, partial differential equations with singular coefficient functions, the 
question how the waves behave in cloaking devices near the surface where 
the material parameters are singular is complicated. Indeed, very different 
kind of behaviors of solutions are possible: In the three dimensional case, 
it is proved in [10] that the transformation optics construction based on 
a blow up map allows cloaking with respect to time-harmonic solutions of 
the Helmholtz equation or Maxwell's equations as long as the object being 
cloaked is passive. In fact, for the Helmholtz equation, the object can be 
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an active source or sink. Moreover, in [10] it is shown that the finite energy 
solutions for the Helmholtz equation in the three dimensional case satisfy 
a hidden boundary condition, namely waves inside the cloaked region sat- 
isfy the Neumann boundary conditions. For Maxwell's equations, the finite 
energy solutions inside the cloaked region need to have vanishing Cauchy 
data i.e., the hidden boundary conditions are over-determined. This leads 
to non-existence of finite solutions for Maxwell's equations with generic in- 
ternal currents [10]. Physically, this non-existence results is related to the 
so-called extraordinary boundary effects on the interface between the cloaked 
and uncloaked regions [39]. 

Another point of view in dealing with the singular anisotropic design for 
cloaking devices is to approximate the ideal cloaking parameters by nonsin- 
gular, or even nonsingular and isotropic, parameters [12, 13, 14, 19, 20, 27], 
which has its advantages in practical fabrication. In the truncation based 
nonsingular approximate cloaking for three dimensional Helmholtz equation 
[14], when it approaches the ideal cloaking, one can obtain above Neumann 
hidden boundary condition for the finite energy solution. In the nonsingular 
and isotropic approximate cloaking, one can obtain different types of Robin 
boundary conditions by varying slightly the way how the approximative cloak 
in constructed, see [12, 13, 14]. 

Similarly, for Maxwell's equations it has been studied how the approximate 
cloak behave on the limit when the approximate cloaks approach the ideal 
one [27]. We note that for Maxwell's equations there are various suggestions 
what kind of limiting cloaks are possible in three dimensions. These sug- 
gestions are based on constructions where additional layers (e.g. perfectly 
conduction layer) is attached inside the cloak [10] or where the ideal cloak 
corresponds to some of the possible self-adjoint extensions of Maxwell's equa- 
tions [36, 37]. In the two dimensional or cylindrical cloaking construction for 
Maxwell's equations, the eigenvalues of permittivity and the permeability of 
the cloaking medium do not only contain eigenvalues approaching to zero 
(as in 3D) but also some of the eigenvalues approach infinity at the cloaking 
interface. Then, the electric flux density D and magnetic flux density B may 
blow up even when there are no sources inside the cloak and an incident plane 
wave is scatters from the cloak, see [11]. However, if a soft-hard (SH)-surface 
is included inside the cloak, the solutions behave well. 

These above examples show how different behavior the solutions may have, 
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in different type of cloaking devices, near the interface between the cloaked 
and uncloaked regions. 

In this paper, we analyze the two dimensional cloaking for Helmholtz equa- 
tion when there are sources present inside the cloaked region. We start with 
the nonsingular approximate cloaking based on the truncation of the singu- 
lar transformation. Taking the limit when the approximate cloaks approach 
the ideal cloak, we show that a non-local boundary condition appears on 
the inner cloak interface. This type of boundary behavior is very different 
from that the solutions have in three dimensional case discussed in [14, 39]. 
The main result is formulated as Theorem 3.2. We note that in the recent 
preprint [31] of Hoai-Minh Nguyen a different type of formulation, based on 
a transmission problem, is given for the non-local boundary condition ap- 
pearing in two-dimensional cloaking. In [31], also cloaking for more general 
second order equations and quantitative convergence properties of the three 
and two dimensional approximative cloaks are analyzed. 

Physically speaking, the non-local boundary condition is possible due to the 
fact that the phase velocity of the waves in the invisibility cloak approaches 
infinity near the interface between the cloaked and uncloaked regions, even 
though the group velocity stays finite, see [7]. We note that as the most 
important experimental implementations of invisibility cloaks [33] have been 
based on cylindrical cloaks, the appearance of such boundary condition could 
also be studied in the present experimental configurations, at least on micro- 
wave frequencies. We also study the eigenvalues, i.e., resonances inside the 
ideal cloak corresponding to the non-local boundary condition. As these 
eigenvalues play an essential role in the study of almost trapped states [13] 
and in the development of the invisible sensors [2, 15], such resonances can 
be used to study analogous constructions in cylindrical geometry. 

The rest of the paper is organized as following. In Section 2, we present 
the basics on transformation optics in the electrostatics setting and apply 
them to the construction of acoustic ideal cloak for the two dimensional 
Helmholtz equation. Section 3 is devoted to the nonsingular approximate 
acoustic cloaking construction and analysis of behaviors of acoustic waves as 
it approaches the ideal cloaking. 
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2 Perfect acoustic cloaking 



2.1 Background: electrostatics 

Our analysis is closely related to the inverse problem for electrostatics, or 
Calderon's conductivity problem [3, 5, 29, 30, 35]. Let Q C M. d be a domain, 
at the boundary of which electrostatic measurements are to be made, and 
denote by cr(x) the anisotropic conductivity within. In the absence of sources, 
an electrostatic potential u satisfies a divergence form equation, 

V-(jVu = (1) 

on Q. To uniquely fix the solution u it is enough to give its value, /, on the 
boundary. In the idealized case, one measures, for all voltage distributions 
u\an = f on the boundary the corresponding current fluxes, u-aVu, where 
v is the exterior unit normal to dfl. Mathematically this amounts to the 
knowledge of the Dirichlet-Neumann (DN) map, A CT . corresponding to cr, 
i.e., the map taking the Dirichlet boundary values of the solution to (1) to 
the corresponding Neumann boundary values, 

K ■ u\ d n h> u-aVu\ dQ . 

If F : Q, — > fl, F = (F 1 , . . . , F d ), is a diffeomorphism with F\gn 
then by making the change of variables y = F(x) and setting u 
we obtain 

V ■ aVv = 0, 
where a = F*a is the push forward of a in F, 

This can be used to show that 

Thus, there is a large (infinite-dimensional) family of conductivities which 
all give rise to the same electrostatic measurements at the boundary. This 
observation is due to Luc Tartar (see [21] for an account.) Calderon's inverse 



Identity, 
v o F' 1 , 

(2) 



5 



problem for anisotropic conductivities is then the question of whether two 
conductivities with the same DN operator must be push-forwards of each 
other. There are a number of positive results in this direction in two dimen- 
sions [4, 22, 25, 34] , but it was shown in [16, 17] in three dimensions and in 
[20] two dimensions that, if one allows singular maps, then in fact there are 
counterexamples, i.e., conductivities that are undetectable to electrostatic 
measurements at the boundary. 

From now on, we will restrict ourselves to the two dimensional case. For 
each R > 0, let Br = {\x\ < R} and E# = {\x\ = R} be the central ball 
and sphere of radius R, resp., in M 2 , and let O = (0, 0, 0) denote the origin. 
To construct an invisibility cloak, for simplicity we use the specific singular 
coordinate transformation F : M 2 \{(3} — > M 2 \i? 1 , given by 

f V, for \y\ > 2, 

* = *W:=| + £orO<M<2. < 4 > 

Letting ao = 1 be the homogeneous isotropic conductivity on I 2 , F then 
defines a conductivity a on M 2 \5! by the formula 

a* k (x) := (F*a y k (x), (5) 

cf. (3). More explicitly, the matrix a = [cr J,fc ] 2 fc=1 is of the form 

a(x) = Jiz^n(x) + — pi— (J-n(x)), K|x|<2, 



where U(x) : 1R 3 — > M 3 is the projection to the radial direction, defined by 

U(x)v = (v~) A, (6) 

V \ X \J fI 

i.e., n(x) is represented by the matrix {x^xx*, cf. [20]. 

One sees that <j(x) is singular at Ei, that is, in interface between the cloaked 
and uncloaked regions, as one of its eigenvalues, namely the one correspond- 
ing to the radial direction, tends to 0, and the other tends to oo as \x\ 
approaching 1 + . We extend a to B\ as an arbitrary smooth, nonsingu- 
lar (bounded from above and below) conductivity there. Let Q = the 
conductivity a is then a cloaking conductivity on Q, as it is indistinguish- 
able from cr , vis-a-vis electrostatic boundary measurements of electrostatic 
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potentials (treated rigorously as bounded, distributional solutions of the de- 
generate elliptic boundary value problem corresponding to a [16, 17]). 

A similar construction based on a blow up map F was proposed in Pendry, 
Schurig and Smith [32] to cloak the region B\ in M 3 from observation by 
electromagnetic waves at a positive frequency; see also Leonhardt [23] for 
a related proposition for the Helmholtz equation in M 2 based on the use of 
several leaves of an Riemannian sufrace. 



2.2 Ideal cloaking for the Helmholtz equation with in- 
terior sources 

We consider the Helmholtz equation, with source term p, of the form 

AV • ffVw + co 2 u = p(x) on f2 (7) 

corresponding cloaking medium with the inverse of the anisotropic mass den- 
sity and the bulk modulus given by 



^■jk 



where (a , A ) corresponds to homogeneous vacuum space and (a a , X a ) are 
arbitrary smooth, nondegenerate medium in cloaked region B\. The push- 
forward of tensor F*o"o is defined by (3) and -F*Ao by 

(F*\ )(x) := [det(DF)X ]oF- 1 (x). 

More specifically, if (cxo, Ao) = (/, 1), one has for 1 < \x\ < 2, 

a {x) = M^n(x) + 1 -pL(/-n(z)) J X(x) 



jk 


for 


\x\ > 2 




[ A 


for 


\x\ : 


(F*a o y k 


for 


1 < \x\ < 2 , 


H 


-F*A 


for 


1 < 




for 


\x\ < 1 




[ A a 


for 


\ x \ '■ 





X 




4(|z 


1-1) 



are both singular at the cloaking surface S := Ei = dB\. 

In the next section, we study the regularized approximate cloaking and obtain 
the behavior of acoustic waves in a singular medium by taking the limit of 
waves propagating in the nonsingular medium. 
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3 Nonsingular approximate cloaking for the 
Helmholtz equation with interior sources 



At the present time, we assume that a and A be homogeneous isotropic in 
side £>!, i.e., (a, A) = (a a 5^ k , A a ) with a a and A a arbitrary positive constants. 

To start, let 1 < R < 2, p = 2(R— 1) and introduce the coordinate transfor- 
mation F R :R 2 \B p ->R 2 \B R , 

( y, for \y\ > 2, 

X : = F ^> = { (l + M) JL for , < \y\ < 2. 

We define the corresponding approximate mass tensor a R and bulk modulus 

\ R as 



a^ k {x) for \x\ > R, , , , _ J X(x) for \x\ > R, 
a a 5 jk for \x\ < R. R ^ X} = \ X a for \x\ < R, ^ ] 

where a^ k and A are as in (8). Note that then a^ k (x) = ((-Fr)* uqY ' (x) 
and X(x) = ((Fr^Ao) (x) for \x\ > R. Observe that, for each R > 1, the 
medium is nonsingular, i.e., is bounded from above and below with, however, 
the lower bound going to 0, and the upper bound going to oo as R —y 1 + . 
Consider the solutions of 

(XrV- ctrV + uj 2 )u r =p in tt ^ 
u R \an —f, 

As ctr and Xr are now non-singular everywhere on Q, we have the standard 
transmission conditions on H R := {x : \x\ — R}, 

Ur\x r + = u R \x R -, 
e r - (TrVu r \x r + = e r - ctrVur\x r -, 

where e r is the radial unit vector and ± indicates when the trace on H R is 
computed as the limit r — y R^ 1 . 

Let Q = .63. Then ur defines two functions such that 



u R (x) 



Vr(Fr\x)), for R < \x\ < 3, 
v R (x), for \x\ < R, 
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and v]i satisfy 

(V 2 + w>+(y) = p(F R (y)) in p < \y\ < 3, 
Vr\8B 3 = f, 



(12) 



and 

(V 2 + k 2 uj 2 )v r (x) = K 2 p(x), in |x| < R. (13) 

where k 2 = (<J a \ a )~~ 1 is a constant. Moreover, if we assume u 2 is not an 
eigenvalue of the transmission problem, then by the transformation law we 
have 

e r -a R Vu R \ m = e r -Vvt\ m . 

This implies that the DN-map A aR: \ R at dQ for the approximate cloaking 
medium (9) is the same as the DN-map at dfl, denoted by A p , of a nearly 
vacuum domain with a small inclusion present in B p . 

Next, using polar coordinates (r, 9), r = \y\, and (r,9), r = \x\, the trans- 
mission conditions (11) on the surface S R yield 



pd r v R (p,9) = KRdpv R (R, 



(14) 



We consider the source term n 2 p(x) where p(x) G C°°(IR 2 ) with suppp C B Ro 
(0 < Rq < 1). It generates a radiating wave w(x) G C 00 ^ 2 ), namely the 
solution of 

(V 2 + «V)iu = K 2 p in M 2 (15) 
satisfying the Sommerfeld radiation condition. Moreover, one can write w as 

oo 

w(r,6) = PnH$(Kur)j n(> , for r > Rq (16) 

n=— oo 

where H$(z) and J\ n \(z) denote the Hankel and Bessel functions, see [8]. 



In Bj^Brq the function v R [x) differs from w by an entire solution to the 
homogeneous equation of (13), and thus for r G (Rq,R) 



(a n J\n\(Kur) +p n H^(nu7))e in0 , 
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with yet undefined a n = a n (K,u; R). Similarly, for p < r < 3, 

oo 

<(r,9) = {cnHQ(ur) + b n J lnl (ur))e md , (17) 

n=— oo 

with as yet unspecified b n = b n (K, u; R) and c n = c n (n, u; R). 
Rewriting the boundary value / on dQ as 

oo 

/(*) = E ^ eme ' 



for a n , b n and c, 



we obtain, together with transmission conditions (14), the following equations 

f n = b n J H (3u) + c n H^(3u), (18) 

) 
I 



a n J\ n \(KuR) + p n Hr)(KuR) = b n J\ n \(up) + c n #, ( 7 (up) , (19) 



KR(nuja n (J\ n \)'(KuR) + KupJHr})' '(kujR)) 

m (20) 
= J| n( ) (cup) + wc^i^) (wp)). 

Solve for a n and c n from (19)-(20) in terms of p n and b n , and use the solutions 
obtained and the equation (18) to solve for b n in terms of /„ and p n . This 
yields 



■(/„ + s n Hr((3u)p ri 



where 



bn = " m (f + « 

J H (3u) + s n H$(3u)) 
Q>n t n b n t n p n 



— {pJ ]n \(KuR)J' n \(up) - K 2 RJ[ rA (KuR)J\ ri \(up)} , 
i- {pH^{up)J' H {up) - p(H^y(up)J H (up)} , 
-t- \k 2 R(H^)'(kujR)J h (kujR) - n 2 RJ[ nl (KuR)H^ 



(21) 



tn = ^~ {k 2 RH{»(u P )(hU)\kuR) - p{H^)\up)H^(KuR)) 
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with D n the common denominator given by 



D n = K 2 RJU^R)H^Aup) - P J H (KuR)(H^y(up). (22) 



3.1 Resonances inside the cloak 

Suppose that the boundary data vanishes, i.e., / = 0. Then by (21), we have 

On _ ,_ , (1)/- - 



J H (3u) + s n H„ (3 



Therefore, one can show 

(t n s n - t n s n )H^(3u) - t n J\ n \{3u) 

J M (3u) + s n H$(3u>) 

_ ^R{H^)'{KuR)h - pH^KuR)^ _ = K v (23) 
pJ\ n \(KujR)l 2 - K 2 RJ' n ^K0jR)h n ' B n 

where 



h = J\ nl (up)H$(3u) - H^(up)J H (3 



|n| V^; ±A \n\ 



(24) 



h = f H (up)H^{3u) - (H^)'(cup)J w (3cu). 

Since for small arguments < x <C 1, 

1 / ' x\ n n = 

U*) n ~ h ~ { ^ (I)"" 1 n > 1, 

\ (|) n n > i ; > W \ m\2 n x- n - x lTX n > 1 



(25) 

when n > 1, where we denote f~g if f — g = o(g) as x — >• 0, and 

i9 n iii~ n ~ 1 (r) — 1 V 
A„ ~ — i -J„(3u;) [ujk 2 R{H^)\kujR) + (/«<;#)] 

— i2 n w" n ~ 1 fra — IV 
5 n i -Jn(3o;) \uK 2 RJ' n (KUjR) + nJ n («;u;#)l /T n , 

7T 
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for n > 1. Now we observe that \a n \ — > oo as R — > 1 + if 



[wk 2 /J( J| n |)'(Ko;i2) + \n\J\ n \(nu)R 



0. 



j?=i 



Note that then 



[uk 2 R(H^)'(kuR) + \ti\H^(kuR 



7^0. 



(27) 



(28) 



i?=i 



This implies that if k is outside a discrete set and if u is such that (27) and 
(28) are satisfied by functions J\ n \ and for some n, then there are sources 
p for which the if 1 -norm of the solution ur goes to infinity in the cloaked 
region (i.e., when resonance happens) as R — > 1 + (i.e., p — > 0). We remark 
that condition (27) implies condition (28) automatically and is equivalent to 
that the function 



V± n (r,6) := J H (Kur)e MnB 
satisfies the boundary value problem 

(A + kV)V = in B u 



[ K rd?V + {-d 2 e f/ 2 V] 



0. 



(29) 



(30) 



F=1H 



Next we consider the frequencies uj for which 



\uik RJ'^(kujR) + \n\J\ n \(nuR)\ 



R=l 



for any n G Z. (31) 



J| n |(3w) ^ 0, 



3.2 Non-local boundary condition with non-resonant 
frequencies 

In the following, we show that when we have in B 2 \Br, 1 < R < 2 the ap- 
proximative cloaking material parameters, then for the non-resonant frequen- 
cies oj, the boundary condition in (30) holds for all solutions when R — > 1 + . 

Lemma 3.1. Assume that in Q = B3 we have the material parameters 
((Tr,Xr). Moreover, suppose u is such that (31) holds. When R > 1 is 
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sufficiently close to 1, then for any source p £ L 2 (Q) supported compactly in 
Bi and for f £ H l l 2 (dVL) the Helmholtz equation (10) has a unique solution 
ur. Moreover, as R —> 1 + , Fourier coefficients of the solution v R := u r \b r 
in the cloaked region, 



"2tt 



the limits lim^i-j. v Rn (r) exists and we have for all n £ Z 

Hm («r^ n (r) + \n\v^ n (r))\ 7=R = 0, (32) 



tim(KrdrV Rn {r) + \n\v R>n (r))\ ¥=1 = 0. (33) 
Proof. We suppose that u is not an eigenvalue of (30). Then 

oo 

««M)= (a n J\n\(^r)+p n H^(^)e me 

n=—oo 
oo 

E 



m0 



A n (R) = K 2 R{H^)\KuR)h - P H^{KwR)k 



where as in (23) 

A fP\ — .-2 or T- 

B n (R) = pJ\ n \(nuR)l2 — K 2 RJl n \(KLuR)li. 

In following, we sometimes denote shortly A n (R) = A n and B n (R) = B n . 
Denote 

*n(r) := ^J|n|(Aewr) + ifft(«wr). 

Apparently $(r, 0) := <J> n (?)e m6 ' satisfies the Helmholtz equation 

(A + kV)$ = in B R . 

To prove the existence of the limits \\m R ^ 1+ v Rn (j r ) and (32), it is sufficient 
to show 

rK^$„(r) + n$ n (r)| F=Jj , ^ as R -> 1 + . (34) 
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Indeed, by 



and lim^i B n (R) ^ (corresponding to the non-resonant case), we have 



(36) 



rKd?$ n (r) +n$ n (r)| F= 

= ^(nh+uph) ({H^)'(KLjR)J\ n \(KuR) - J' lnl (KUjR)H^(KLoR^ 

tin ^ 

where l\ and Z 2 are given by (24). 
For n > 1, as p — > + , 

nh + wp/ 2 =^P f J^-i^p^^^a;) - ^^(w/o) J| ft |(3o; 
=0{p~ n+2 ). 

Combining (35), (36) and (26), one has 

TKd& n (j) + \n\$ n (?)\ ¥=R = 0(p 2 ) as R 1 + (i.e. p -)• + ), (37) 

which proves (34). 

For n = 0, from (25), one has 

Ao = - ^f^( y Hi ) 1) )'(KuR)J (3u) In (^) + k 2 R(H^)'(kujR)H^ (3w) 

7T V 2 / 

2z 

+ —H£\kuR)J (3uj) + 0{p), 



TXUO 

-2 ; 



B =^-^^(«wi2)J (3o;)ln - k 2 RJ' (kuR)H^\3 

it V 2 / 

2z 

+ — J (kujR)J (3cu) + 0(p). 

7VUJ 

Therefore, 

5^ (i2) = ^-kuJ'o{kuR) + Ku(H^y(K0jR) (38) 

-DO 

has denominator i?o and numerator 

«a;[A Jo(«a;E) + Bo(H^Y(kujR)] = —W n (KuR) + £>(p), 

7T 
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where W n (x) = Hq 1 \x)Jq(x) — (Hq)'(x)Jo(x). This implies 

drMR) W^KUK) 

V ; kRJ^kuR) In (f) ' 

i.e., the boundary condition (34) is satisfied for n = and moreover, 

rKdrMr) \ r=R = O \ — ^ ] as i2 1 + . (39) 

This proves (32). The equation (33) follows similarly by evaluating (35) and 
(38) at f = 1 instead of r = R. □ 

Now we are ready to prove our main result of the limit of the waves ur of the 
physical approximate cloaking medium as R — > 1 + . We recall that E = dB\. 

Theorem 3.2. Let u be such that (31) is satisfied. Assume that ur is the 
solution of (10) with f = andp G C^°(Br ) with R < 1. Then as R — >■ l + ; 
u# converges uniformly in compact subsets of B 3 \H to the limit Ui satisfying 

(V 2 + K 2 bj 2 )u x = k 2 P in Bx, (40) 

KdrU! + (-^) 1/2 Ul \ dBi =0, (41) 

and 

w iL\^ = - ( 4 2) 

We note that solutions of (40)- (42) with / ^ and p = have been analyzed 
in [11]. 

Proo/. Let R < Ri < R. Recall that solution w E C°°(IR 2 ) of (15) is the 
radiating solution produced by source K 2 p in M 2 and it has the expansion 
(16) for r > R\. Consider the Fourier coefficients 

/•2tt 

Wn (r) = / e- ine w(r,6)d9 
Jo 

and denote P n (x) := —(V 2 + k 2 — n 2 )w n (\x\). As w G C 00 ^ 2 ), we see using 
integration by parts that 

\\Pn(\x\)\\ L i(B 3 ) < C M (1 + \n\y M for arbitrary M > 0. (43) 
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We consider the Fourier coefficients 



2tt f2ir 



t£ B (r)=/ e-"»vJi(r,e)dB, «J B (r) = / e-» %+(r,0)eZ0. 
Jo Jo 

They satisfy the following problem 

(-V 2 -KV + n 2 )^ n (|x|) = P„(|x|) forO<|x|<i?, (44) 
(-V 2 - a; 2 + n>+ n (M) = for p < |y| < 2, (45) 
4nU = 0, (46) 

U fl,nlaB- = Wfl,Jsfl+» ^(^J'G^DU" = PK,J'(M)leB+- ( 47 ) 

By the transmission conditions (47), we see for V Rn (x) = v Rn (\x\) 

R f p 



drV^ n ds x = - / -^d r v+ n v+ n ds y = / -drV+y^dSy. 

dB R P JdB p KIt JdB p K 

Thus, using integration by parts, we obtain 

h -=\ [ PnV^dy 
Jb r 

= [ ("V 2 - «V + n 2 )V^ n V^ n dx + - f _(-V 2 - u 2 + n^V+J^dy 
Jb r k Jb 2 \b 



dB R \JdB 2 JdB i 



drV^ n dS x - \ I - I ] -d r V+ n V+ n dS y 



[ (\^ n \ 2 + (-K 2 u 2 + n 2 )\V^ n \ 2 )dx+^ I _(\VV+J 2 + (-u 2 + n 2 )\V+ n \ 2 )dy, 

J Br K JB 2 \B p 

and then 



h=\ (|V^- n | 2 + (-/,V + n 2 )|V R - n | 2 )rfx + - / _(|Vy+J 2 + (-o; 2 + n 2 )|V R + ri | 2 )cZ 2/ 

J Br K JB 2 \B p 

>[ (-KW + n 2 )\V^\ 2 dx + - [ _(-u 2 + n 2 )\V+ n \ 2 dy. 

J Br K JB 2 \B p 

For \n\ > N > max{K 2 u 2 ,w 2 }, the above and Ii < \\V Rn \\ L 2\\P n \\ L 2 imply 
first that 

(\\ V R,n\\LHB R ) + II^JIz, 2 (B 2 \B7)) < C N \\ Pn || L 2 (B 3 ) ( 48 ) 
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and second that 

(II^JIj^CBh) + \\ V R,n\\m(B 2 \B;)) < C' No \\P n \\ L 2 {B3) (49) 
where Cn and C' N are independent of R and n. 

By Lemma 3.1, for each n G Z and r G [0, 1) and r G (0, 2) there exists limits 
v n (r) = Jim (r) > ( r ) = Jim v£ n (r) (50) 

and we denote V^(x) = v^(\x\). Let now < r± < R\ < 1. Then by (49) 
the restrictions ^Jo-i,/?!), R> R\ are uniformly bounded in H 1 (r 1 , Ri). By 
Sobolev embedding theorem, the set {^fl re |[ri,Ki]; -Ro < R < 1} is relatively 
compact in C([ri,Ri]) C H s (ri, R±), 1/2 < s < 1. Thus any sequence 
{ v Rj n\[ri,Ri])j^=i wr th Rj — > 1 has a subsequence converging in C([ri,Ri]) 
which limit has to coincide with v~ by (50). Thus v Bn have to converge to 
v~ in C([ri,Ri\) and hence V Bn have to converge to V~ in C{B Bil \ B ri ) 
as R —t 1. Similarly, for all p\ > we see using (49) that V Rn have to 
converge to in C(B 2 \ B P1 ) as -R — > 1. Now, as the Sobolev norm u 
1 1 u 1 1 if 1 (Br \B r ) * s a l° wer semi-continuous function in L 2 (B Bl \ B n ) and (49) 
holds, we see that \\V~\\ H i^ BR ^ Bri) < C' No \\P n \\ L 2 (Ba) for all < r x < R x < 1. 
Hence, by (50) we see using e.g. monotone convergence theorem and [18] that 

IlKTll^HBi) = IIK~IIhi(Bi\o) = ri Jjm^ 1 llK"llHi(B Hl \flr 1 ) - ^kll^lUw-t 51 ) 

By (43) we see that Wi(r, 6 1 ) = Xlnez (r)e m0 is a well defined function in 
H 1 (B 1 ) satisfying (40). By (49), 

ll^nllc(B fll \B n ) = II^JIcCK-Ri]) < ^N 0l n,Ri \\ v R,n\\H x (rx,Ri) ^ ^Wo.n.Jii l|-fn|U 2 (B 3 )' 

where Cjy n fll does not depend on i? or n > iVo- Thus, using (43) and the 
convergence of V^~ n to in C(B Rl \ B n ) we see that u R converge to u\ 
uniformly in compact subsets of Bi \ {0}. Using equation (44), we see the 
uniform convergence in a neighborhood of zero, too. Similarly, we see the 
uniform convergence in compact subsets of B 3 \ B\. Moreover, by [26] and 
(40) the boundary values K,d r Ui\g Bl G H~ l l 2 (dBi) and u 1 \g Bl G ^^(dBi) 
are well defined. Now, by (44) and [26] we have \\nd r V~\[ H -yzf dB ^ < C 2 (l + 
n ' 2 )\\Vn\\H 1 (B 1 \B Ro ) where C2 > is independent of n. Thus using (43), (33), 
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and (51) we see that the boundary value nd r Ui + (— dgY^Ui, vanishes. Hence, 
the boundary condition (41) is satisfied. 

Equation (42) for Ui\ Bs ^g^ follows using (43), (49) and the fact that by (21), 
(22) for any fixed n G Z we have in (17) 

c n (p) = 0(pW), b n (p) = O(pM) asp^0+, |n|>l, 
c b (p) = 0((logp)- 1 ), 6 o (p) = 0((logp)" 1 ) asp^0+. 

□ 

We conclude our discussion by remarking that our analysis also explains 
the limit of approximate electromagnetic cloaks in the cylindrical case with 
TE/TM polarized incoming waves, as the solutions of Maxwell's equations 
in this case satisfy the Helmholtz equation. 
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